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We report inelastic neutron scattering measurements of the spin dynamics in the layered hexagonal 
magnet 2//-AgNi02, which has stacked triangular layers of antiferromagnetically-coupled Ni^+ spins 
(S'^l) ordered in a coUinear alternating stripe pattern. We observe a broad band of magnetic 
excitations above a small gap of 1.8 meV and extending up to 7.5 meV, indicating strongly dispersive 
excitations. The measured dispersions of the boundaries of the powder-averaged spectrum can be 
quantitatively explained by a linear spin-wave dispersion for triangular layers with antiferromagnetic 
nearest- and weak next-nearest neighbor couplings, a strong easy-axis anisotropy and additional 
weak inter-layer couplings. The resulting dispersion relation has global minima not at magnetic 
Bragg wavevectors but at symmetry-related soft points and we attribute this anomalous feature 
to the strong competition between the easy-axis anisotropy and the frustrated antiferromagnetic 
couplings. We have also calculated the quantum corrections to the dispersion relation to order 
1/S in spin- wave theory by extending the work of Chubukov and Jolicoeur [Phys. Rev. B 46, 
11137 (1992)] and find that the presence of easy-axis anisotropy significantly reduces the quantum 
renormalizations predicted for the isotropic model. 



PACS numbers: 75.25.-(-z, 75.10.Jm 
I. INTRODUCTION 



Quantum antiferromagnets on triangular lattices pro- 
vide model systems for investigating the effects of quan- 
tum fluctuations and geometric frustration. Zero point 
fluctuations are expected to be strongly enhanced for low 
spin and frustrated lattice geometry and may stabilize 
non-classical ordered or spin liquid phases, or unconven- 
tional spin dynamics. The triangular Heisenberg anti- 
ferromagnet with first and second neighbor couplings Ji 
and J2 [see Fig. [1] shows strong frustration effects and 
a macroscopically-degenerate classical ground state for 
1/8 < J2/J1 < 1. This was initially proposed theoreti- 
cally as a candidate for a chiral spin liquid state. ^ How- 
ever, perturbative expansions using a spin-wave basis^'^"^ 
predicted that in this range of couplings quantum fluctu- 
ations do not stabilize a spin liquid state but instead lift 
the classical degeneracy to select a collinear stripe order 
pattern [see Fig [T]. Moreover, large quantum renormal- 
izations of the spin- wave dispersion relations compared to 
classical have been predicted, but not probed experimen- 
tally in the absence of a good experimental model sys- 
tem. The recent observation of a collinear stripe-ordered 
ground state in the layered hexagonal antiferromagnet 
2i?-AgNi02 [ref. ^ suggests a potential realization of 
the experimentally-unexplored frustrated Ji — J2 model 
in the collinear phase and here we present flrst inelas- 
tic neutron scattering measurements which give informa- 



tion about the spin gap and dispersion relations. We 
find that the data can be parameterized by a spin-wave 
dispersion relation for a J1-J2 triangular antiferromag- 
net with strong easy-axis anisotropy and weak inter-layer 
couplings. 

Delafossite materials of the type XNi02 have been gen- 
erally thought of as possible two-dimensional frustrated 
magnets. A network of edge-sharing NiOe octahedra lead 
to a triangular lattice arrangement of Ni ions in planes 
spaced by layers of ions. There have been a vari- 
ety of studies investigating the properties of compounds 
where X=Li, Na. LiNi02 has shown no long-range mag- 
netic order but experiments are hindered by the diffi- 
culty in preparing stoichiometric samplesi^ From elastic 
and inelastic neutron scattering measurements, NaNi02 
was found to be a spin-1/2 system with in-plane (un- 
frustrated) ferromagnetic interactions and weak antifer- 
romagnetic inter-layer couplingsj^ 

We have recently started exploring the delafossite 
AgNi02, which shows frustrated antiferromagnetic in- 
plane interactions^'^. Detailed structural studies have 
been performed on the hexagonal polytype of AgNi02, 
so called 2//-AgNi02 with two Ni02 layers per unit 
cell [see Fig. [^b)], as opposed to the earher-synthesized 
rhombohedral 3i? polytype^ with a 3-layer stacking se- 
quence along c-axis. High-resolution neutron diffraction 
measurements of 2i7-AgNi02 have observed a structural 
transition upon cooling below 365 K, which leads to a 
periodic arrangement of expanded and contracted NiOg 
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FIG. 1: (Color online) Phases of the Heisenberg triangular 
lattice antiferromagnet (120° coplanar order, collinear stripes 
and incommensurate spiral) as a function of the ratio between 
the first and second neighbor couplings, Ji and J2 respec- 
tively. In the collinear phase the common spin direction is 
spontaneously chosen and can point anywhere, it is shown 
here in-plane for ease of visualization. 
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FIG. 2: (Color online) a) Triangular lattice formed by Nil 
sites (dark brown spheres) in a hexagonal Ni02 layer in 2H- 
AgNi02; thick arrows indicate the paths for the in-plane nn 
and nnn exchanges Ji and J2. Dashed diamond shows the 
crystallographic unit cell. Gray spheres are Ni2 and Ni3 sites, 
assumed non-magnetic, b) 3D view of the crystal and mag- 
netic structure. There are two symmetry-equivalent Ni02 lay- 
ers per unit cell. Layer stacking is such that Nil ions sit above 
the centre of a Nil triangle in the layer below. We consider 
two natural candidates for inter-layer couplings: Jz between 
Nil ions in adjacent layers (three neighbours above and three 
below) and between Nil ions at two layers apart. 



octahedra.^^ This was proposed^ to be a consequence of 
spontaneous charge order on the Ni sites driven by the 
need to lift a two-fold orbital degeneracy. This leads to 
a strongly-magnetic Nil site (spin-1 Ni^+, l/'S^'^ of sites) 
arranged in an ideal triangular lattice [see Fig.[2k)], with 
the remaining 2/3''^'' of Ni sites being Ni'^-^+ with an itin- 
erant character. 

Magnetic order occurs below T/v = 19.7(3) Ki^ when 
the magnetic moments of Nil sites order in a collinear 
pattern of alternating ferromagnetic rows (stripes) in the 
triangular layers with spins pointing along the c-axis, as 
shown in Fig. [3^). A large ordered moment is found on 
the Nil sites, 1.552(7) at 4 K, and no significant or- 
dered moment could be detected on the remaining Ni2 
and Ni3 sites, proposed by band-structure calculations'^ 
to be strongly itinerant and possibly ordered but with 
only a very small (< O.l/ie) moment. So to a first approx- 
imation the coherent spin dynamics in the magnetically- 
ordered phase is expected to be dominated by the S = 1 
Nil spins arranged on an ideal triangular lattice with 
antiferromagnetic couplings. 

The experimentally observed collinear stripe order in 
Fig. [3^) is rather unusual for a triangular antiferromag- 
net (TAFM), but has been proposed to occur for the 
classical Isingmodel for finite second-neighbor couplings 
J2 > [Ref. [il|. For the Heisenberg TAFM the collinear 
stripe order has also been proposed to occur as a ground 
state for moderate J2 in the range 1/8<J2/Ji<l, Ref. 
1^. In this range the classical ground state is macroscopi- 
cally degenerate and there are many other non-coUinear 
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FIG. 3: (Color online) a) Magnetic structure in a triangular 
layer showing the alternating stripe order of the Nil spins (± 
symbols indicate the projection of the ordered spin moments 
along the c-axis). x, y labels indicate a natural orthogonal 
coordinate system with x along and y transverse to stripes. 
Solid diamond indicates the magnetic unit cell, doubled along 
the crystallographic a-axis compared to the chemical unit cell 
(light gray shaded area). The depicted magnetic structure 
has stripes running along the b-axis and ordering wavevector 
fc = (1/2,0,0); equivalent domains are obtained by rotating 
this pattern by ±60° around the c-axis. Further neighbor 
couplings are required to stabilize this structure with respect 
to the ferrimagnetic honeycomb up-up-down pattern shown 
in b). 



states degenerate with the two-sublattice collinear stripe 
order shown in Fig. [3^), but zero-point quantum fluc- 
tuations are predicted to select the latter through the 
"order by disorder" mechanism. The large classical de- 
generacy is manifested also in the linear spin-wave dis- 
persion which has many soft points with zero energy, and 
proper inclusion of quantum corrections leads to renor- 
malizations of the semiclassical dispersion relation and a 
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gapping of the non-physical zero modes.-3. Those effects 
have remained experimentally largely unexplored. The 
recent observation^ of a collinear stripe ordered phase in 
the hexagonal magnet AgNi02 suggested a possible ex- 
perimental realization of the S — 1 TAFM in the range 
of moderate frustration and motivated us to measure its 
spin dynamics in some detail. We note that the physics 
we are exploring here appears to be different to that in 
the recently explored S = 5/2 TAFM CuFe02 which has 
a related, but different collinear structure of alternat- 
ing "double stripes" (ttii) stabilized by strong three- 
dimensional and third-neighbour in-plane couplings In 
2if-AgNi02 the magnetic structure is a (simpler) single- 
stripe pattern t| with only two in-plane sublattices pos- 
sibly stabilized by fewer exchanges (only first and second 
in-plane neighbours) and easy-axis anisotropy. 

The plan of this paper is as follows. The following 
Sec. ini gives details of the inelastic neutron scattering 
experiments to probe the powder-averaged spin dynam- 
ics and the results are presented in Sec. IIIII In the fol- 
lowing Sec. IIV Al the data is parameterized in terms of 
an empirical sinusoidal dispersion model with minima at 
magnetic Bragg wavevectors and different zone bound- 
ary energies along the three orthogonal directions in the 
Brillouin zone; this parameterization gives a gapped and 
predominantly two-dimensional dispersion relation. In 
the following Sec. IIVBI the data is compared with lin- 
ear spin-wave theory for a microscopic spin Hamiltonian 
for the localized Ni spins that includes both first- and 
second-neighbour antifcrromagnetic couplings in the tri- 
angular layers, easy-axis anisotropy modelled by a single- 
ion term, and different models for inter-layer couplings 
consistent with the crystal structure; a minimal spin 
Hamiltonian is proposed. In the following Sec. IIV Cl we 
consider how magnon interactions included at order 1/5 
in spin-wave theory renormalize the dispersion relation 
and find slightly renormalized values for the proposed 
spin interactions. Finally, the results for the spin dy- 
namics are summarized and discussed in the concluding 
Sec.|Vl Technical details of the spin wave calculations for 
coupled easy-axis triangular layers are given in Appendix 
1X1 Details of the derivation of the 1/S quantum correc- 
tions to the dispersion relation for an easy-axis triangular 
antifcrromagnet are given in Appendix [Bl 



II. EXPERIMENTAL DETAILS 

Powder samples of the hexagonal polytype of AgNi02 
were prepared from Ag20 and Ni(0II)2 using high oxy- 
gen pressures (130 MPa) as described in Ref. M- The 
samples used for the inelastic neutron experiments were 
part of the same batch used in previous neutron diffrac- 
tion measurements^^ performed on about half of the total 
sample quantity, which revealed a high-purity hexagonal 
phase (< 1% admixture of the rhombohedral polytype). 
The crystal structure of 2_ff-AgNi02 is shown in Fig.[2h) 
and is hexagonal with space group P6322 (no. 182) and 



lattice parameters a = 5.0908(1)A and c = 12.2498(1)A 
[Ref. [H]. Ni ions are located inside hexagonal (ab) layers 
and there are three distinct crystallographic sites: Nil 
which sit inside slightly expanded NiOg octahedra ar- 
ranged in a periodic triangular lattice of spacing a, sur- 
rounded by a honeycombe of contracted NiOg octahedra 
which contain the Ni2 and Ni3 sites, see Fig. [5^). 

The powder-averaged magnetic excitation spectrum 
was probed using two direct-geometry time-of-flight neu- 
tron spectrometers: MARI, at the ISIS Facility in the 
UK, and IN6, at the Institute Laue-Langevin in France. 
Measurements on MARI showed that the full dynamic 
range of the spin excitations extended only up to 7.5 meV 
so could be accessed with incident neutrons of energy 
Ei = 18 meV, which gave an energy resolution of 0.61(1) 
meV (FWHM) on the elastic line. Higher-resolution mea- 
surements to probe the low-energy part of the spectrum 
were made using IN6 operated with incident neutrons 
of energy Ei = 3.86 meV, which gave a measured en- 
ergy resolution of 0.142(1) meV (FWHM) on the elastic 
line. The sample was cooled using either a closed-cycle 
refrigerator (MARI, base temperature 4.7 K) or orange 
cryostat (IN6, base temperature 1.8 K). Measurements 
were made at the lowest temperatures in the magneti- 
cally ordered phase, near the ordering transition and at 
high temperatures in the paramagnetic phase. 

For the MARI experiment 23 g of powder were placed 
in an annular-shaped sachet to minimize absorption. For 
the low-energy IN6 measurements a similar powder quan- 
tity was placed in a plate-shaped container angled at 
45° to the incident beam and the data was corrected 
for neutron absorption effects using a numerical calcula- 
tion for plate-shaped samples. The scattering intensities 
from both instruments have been converted into absolute 
units of S{Q,U!) of mbarns/meV/sr/Nil by normalizing 
the raw counts to the sample mass and to the measured 
scattering intensities from a Vanadium standard. 



III. MEASUREMENTS AND RESULTS 

An overview of the measured inelastic neutron scat- 
tering data is shown in Fig. |4l Below the Neel ordering 
temperature T/v=19.7(3)K a strong band of scattering is 
clearly observed at low energies, E, and low wavevectors, 
Q = \Q\, [Fig. 3^) and c)], attributed to magnetic excita- 
tions. The intensity of the signal decreases with increas- 
ing wavevector Q, confirming its magnetic character. The 
intensity starts above a small gap of 1.8 meV and ex- 
tends up to 7.5 meV, above which the signal decreases to 
background level. This energy scale for the magnetic ex- 
citations spectrum is consistent with susceptibility mea- 
surements which observed a Curie- Weiss temperature of 
= —9.2 meV— . Energy scans through the data at 
constant wavevectors showing the low energy gap and the 
full extent of the scattering are shown in Figs. [SJ'V and E 
(solid points). The dashed line in those plots indicates 
the estimated non-magnetic background, obtained by a 
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FIG. 4; (Color online) Powder-averaged magnetic excitation spectrum in 2_ff-AgNi02 in the magnetically-ordered phase (top 
panels a) and c)) showing a band of spin excitations above a small gap, b) in the paramagnetic phase at high temperature 
showing an overdamped signal and d) slightly above the magnetic ordering temperature showing a filling-in of the spin gap 
[compare with c)]. Data is raw counts collected after 35 hours (a) and 4.5 hours (c) and normalized to absolute units as 
described in the text. The streak of intensity in panel c) near is non-magnetic and attributed to an acoustic 

phonon dispersing out of the (002) structural Bragg peak. Middle panels show the intensity distribution for the empirical 
dispersion model in eqs. ([1} and (|3}. The calculations include convolution with the instrumental resolution, the magnetic form 
factor, polarization factor and an estimate of the non-magnetic background to be directly compared with the data in the panels 
immediately above. Thick bold arrows labelled A-E indicate the location of scans plotted in Fig. [5^-E and the grey dotted 
lines in all plots indicate the low-Q edge of the measured region. 
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smooth interpolation of the signal observed at energies 
below and above the magnetic signal. The energy gap to 
magnetic excitations fills up upon heating through the 
Neel temperature [compare Fig. St) with d) and scans 
shown in Fig. [6] . The structured band of magnetic scat- 
tering is replaced as expected by an overdamped signal in 
the paramagnetic phase at high temperatures, compare 
Fig. 2^) with b) and various scans in Fig. [Sfl. 

At the lowest temperatures the boundaries of the mag- 
netic scattering show considerable structure which can be 
used to impose constraints on the underlying dispersion 
relation. In particular information is contained in the 
gap, extent of bandwidth, slope of low-Q dispersion up to 
the first minimum gap [see Fig. [4^)] and the low-energy 
dispersion (~ 10% of bandwidth) between subsequent 
minimum gap wavevectors [see Fig. EJ;)]. The powder 
data is a spherical average of the dispersions along all 
directions in reciprocal space weighted by the neutron 
structure factor, so the low-energy boundary of the pow- 
der data corresponds to the dispersion along some direc- 
tion in reciprocal space which has the minimum energy 
at a given Q. In general the spin- wave dispersion has 
global minima at magnetic Bragg wavevectors and here 
strong low-energy scattering is expected. This is con- 
sistent with the data in Fig. \^) showing clear lobes of 
magnetic scattering intensity coming down in energy near 
wavevectors 0.7 and 0.9 where the first two magnetic 
Bragg peaks (1/2,0,0) and (1/2,0,1) occur; a constant en- 
ergy scan near the gap minimum is shown in Fig. [Sp. 
Inelastic intensity occurs near the same wavevectors as 
magnetic Bragg peaks in the diffraction pattern shown 
in Fig. [5p bottom panel [data from the GEM diffrac- 
tometer at ISIS]. The low-energy boundary of the scat- 
tering is clearly modulated as a function of wavevector 
and the onset energy varies as a function of wavevector 
Q as illustrated by comparing energy scans in Fig. [SJ3 
and C. 



IV. ANALYSIS 

The modulations observed in the lower boundary of 
the magnetic scattering indicate strongly dispersive exci- 
tations, however in principle it is not possible to extract 
precise dispersion relations from the powder data, which 
represents an average along all directions in reciprocal 
space. Nevertheless, if the dispersion is highly asymmet- 
ric along different directions and if in-plane and out-of- 
plane lattice parameters are sufficiently different then the 
lower-boundary of the scattering can be identified over 
certain Q-ranges with the dispersion along specific direc- 
tions in reciprocal space, so one may impose certain con- 
straints on the dispersion model. Following this approach 
we have found two models for the dispersion relation that 
could both explain the observed wavevector-dependence 
of the boundaries of the powder-averaged spectrum and 
present them in the following. 



A. Parameterization by an empirical sinusoidal 
dispersion model 

We first consider a parameterization of the data by an 
empirical sinusoidal dispersion model with a finite gap 
A above the magnetic Bragg peaks and different zone 
boundary energies along the three orthogonal directions 
in the Brillouin zone: Ei along the magnetic stripe direc- 
tion, E2 transverse to stripes in the plane, and in the 
inter-layer direction. In detail we consider the following 
empirical form for the dispersion relation ujq 

{uoQf= A2 + (^2_^2) sin2(g,a/2) + 
+ (£;| - A2) sin2(QyaV3/2) + 

+ (El - A2) sin2(Q,c/2) (1) 

where Qx, Qy, Qz sue components (in A^-'^) of the 
wavevector transfer Q in an orthogonal coordinate sys- 
tem with X along the magnetic stripes, y transverse to 
stripes in plane, z normal to the planes. The transfor- 
mation to wavevector components (h, k, I) in units of the 
reciprocal lattice of the hexagonal unit cell is given by 

Qx = fcj 
a 

Qy = ^(2/i + fc), 
av o 

2-77 

Qz = —I. (2) 

c 

Fig. [7] shows a plot of this dispersion relation along 
various symmetry directions in the Brillouin zone and 
Fig. [S}d) shows a 2D contour map of the dispersion. 
The minimum gap A is reached at the origin and at 
all magnetic Bragg peak positions {H, K,L) zt k with 
k — (1/2,0,0) the ordering wavevector and H, K and 
L integers. In order to compare this dispersion relation 
with the data one needs also a model for the neutron 
scattering cross-section. We assume the following simple 
form for the one-magnon cross-section 

^--(Q,,^) = syy{Q,uj) 

= cl-{uk)^,J—^G[u^u:Q) (3) 

Z Lu 

Here C is an overall scale factor and (1— 7q) is a geometric 
factor to concentrate the intensity near the antiferromag- 
netic Bragg peak wavevectors (i/±l/2, L) with iJ, K, 
L integers, 7q — cos((5^a/2) cos((3j^a\/3/2). The 1/uj 
factor is a characteristic energy-dependence of the spin- 
wave structure factor for antiferromagnets and (wfc)jjja,x 
the maximum magnon energy for all wavevectors fc, and 
is introduced here for dimensionality purposes. 5' = 1 is 
the spin quantum number. Eq. ([3]) includes the fact that 
spin fluctuations occur in the plane transverse to the or- 
dered spin direction in a coUinear antiferromagnet, so x 
and y refer to two orthogonal directions in the hexagonal 
ah plane and z is along c. This form is a generic struc- 
ture factor for many simple spin- wave models of coUinear 
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FIG. 5: (Color online) Scans through the powder data along special directions labelled A-E in Fig. ID Scan D shows that there 
is structure in the low-energy inelastic signal for wavevectors where Bragg peaks occur (lower panel in D shows diffraction data 
indicating magnetic Bragg peaks). The energy scan E integrates over a broad Q-range and shows the full bandwidth of the 
magnetic excitations: low temperature data (solid points) is replaced in the paramagnetic phase at high temperatures by an 
overdamped signal (open triangles, solid line is guide to the eye). In panels A-D the thick solid lines show calculations for the 
sinusoidal dispersion model in eqs. ^ and (|3]) with dispersion relations plotted in Fig. [T] The green dashed lines indicate the 
estimated non-magnetic background. 



antifcrromagnets. Here G{u} — loq) is a gaussian function 
which models the instrumental resolution. 

A spherical average of the intensity distribution of eq. 
([3]) including the neutron polarization factor, magnetic 
form factor for Ni^+ ions [as in eq. (|A5p ] and an esti- 
mated additive non-magnetic background (including the 
incoherent elastic line) was compared to the data and 



the overall comparison is shown in the middle panels of 
Fig. m The overall bandwidth and the observed struc- 
ture in the data at low energies is well reproduced by 
this model for gap A = 1.78(5) meV, £'i=5.73(4) meV, 
E2 = 4.15(6) meV and = 2.33(5) meV in eq. The 
maximum magnon energy is (yik)^^^ = 7.0 meV and an 
overall scale factor C ~ 0.35 in eq. ([3]) gives intensities 



Energj' (meV) 

FIG. 6: (Color online) Energy scan at the minimum gap 
wavevector showing a filling of the gap upon heating above the 
magnetic ordering temperature. Solid lines are guides to the 
eye and the dashed line indicates the estimated non-magnetic 
background. 




FIG. 7: Sinusoidal dispersion relation in eq. HJ plotted along 
symmetry directions in the Brillouin zone in Fig. [SJj) (dashed 
lines). Dispersion parameters are A = 1.78 meV, i5i=5.73 
meV, E2 = 4.15 meV and E^, = 2.33 meV. 



comparable to data. 

To analyze the agreement quantitatively we first show 
in Fig. IHK) how the lower boundary of the powder- 
averaged magnetic signal can be identified over certain 
wavevector ranges with the dispersion along particular 
directions in reciprocal space. In particular the observed 
sharp low-Q dispersion up to the first gap minimum near 
Qm = 0.7 in Fig. 0^) is attributed to the disper- 
sion between the origin r(OOO) and the closest magnetic 
Bragg peak at M(l/2,0,0); scan A along the direction in- 
dicated by the bold arrow in Fig. [4^) is therefore very 
sensitive (at low energies) to the FM dispersion and a 
fit to this scan shown in Fig. [5j^ gives E2 = 4.15 ± 0.06 
meV, the zone boundary along the direction transverse 
to stripes in plane. To explain the observed full extent of 
the bandwidth up to 7.5 meV another dispersion, larger 
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FIG. 8: (Color online) a) Reciprocal basal plane of 
2ff-AgNi02 showing the Brillouin zone edges (thick line 
hexagons), zone centers (solid points), and locations (stars) 
of magnetic Bragg peaks from the ordered domain shown in 
Fig. |3^) with stripes along the 6-axis. Labels F, M, K indi- 
cate special symmetry points in the hexagonal Brillouin zone 
(M' and M" are Bragg peak positions for equivalent ordered 
domains with stripes along a-\-b and a, respectively). Dashed 
gray lines shows directions along which the dispersion is plot- 
ted in Fig. [TD] Arrows labelled x, y indicate the natural or- 
thogonal coordinate system used for the empirical dispersion 
model in eq. ((l|. b,c) Color contour maps of the dispersion 
in [ah) plane in the two models consistent with the data: b) 
the empirical parameterization with a sinusoidal dispersion 
plotted in Fig. [7] and c) the linear spin wave model for Ji — J2 
triangular layers with easy axis and inter-layer couplings in 
Fig. llOU Kdashed line). Note that in model b) the minimum 
gap occurs at the magnetic Bragg wavevector M, whereas in 
model c) the minimum gap occurs at the soft point M'. 



in magnitude is required and this is assigned to the direc- 
tion along the stripes, Ei, since the inter-layer dispersion 
bandwidth i?2 — A is assumed to be small. Quantitative 
fits give El = 5.73±0.04 meV. The energy scan in Fig.[5j3 
is sensitive to the gap above the magnetic Bragg wavevec- 
tor, obtained as A = 1.78(5) meV. The energy scan in 
Fig. [5p is at a wavevector corresponding to the midpoint 
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FIG. 9: (Color online) a) Dispersion relations for the empirical 
sinusoidal model in eq. ([1} along various directions involved 
in defining the low-energy edge of the powder-averaged spec- 
trum, b) Corresponding powder-averaged spectrum. Stars 
in both panels are the experimentally-determined low-energy 
edge of the magnetic scattering. Thick arrows indicate loca- 
tion of scans B,C, D in Fig. [5] 



between two subsequent Bragg wavevectors along the c- 
axis, (1/2,0,1) and (1/2,0,2) [see Fig. [SJi)], so this scan 
is sensitive to the inter-layer zone boundary and fits give 
Ez = 2.33 ± 0.05 meV, thus obtaining the last parameter 
of the empirical dispersion relation. 

Fig. [5] shows in more detail the agreement between the 
experimentally-determined onset of the magnetic scatter- 
ing (stars) and the empirical model (various solid lines) 
in panel a) and the calculated powder-averaged spectrum 
in panel b). The onset points were determined from scans 
through the measured powder data at fixed wavevector 
Q and locating the energy where the magnetic intensity 
was above a minimum threshold value. The model repro- 
duces well the shape of the lower edge of the scattering. 
Fits to specific scans shown in Fig. [5l\-D give good ac- 
count of the onset and upper boundary of the scattering. 
The form of the intensity distribution within the excita- 
tion band is not replicated in detail, as seen for example 
in Fig. [5]A_ where the simulation gives a more pronounced 
trough in the intensity in the middle of the band. This 
is almost certainly due to the over simplification of the 
model. 



B. Parameterization in terms of an easy-axis 
Heisenberg model on a triangular lattice 

The empirical sinusoidal dispersion model considered 
in the previous section is the simplest dispersion model 
that can be used to parameterize and explain quan- 
titatively the observed dispersion of the boundaries of 
the powder-averaged magnetic inelastic scattering. This 
analysis suggested that the dispersion relation is gapped 
and predominantly two-dimensional. Here we compare 
the data quantitatively with predictions of linear spin 
wave theory for a microscopic model of coupled trian- 
gular layers where we include antiferromagnetic in-plane 
couplings for the first and second-neighbour exchanges 
Ji and J2, an easy-axis anisotropy and allow for both 
first- and second-layer couplings Jz and J^, as shown in 
Fig. Uti-ib) . Specifically, we consider the spin Hamiltonian 

T^o = X^NN Ji^i ' + Snnn J2Si ■ Si 

J2i i^i) + X^IntcrlaycrNN '^^Si ■ Sk 
+ X^IntcrlaycrNNN "^'z^i ' ^ri, 



(4) 



where the spin is S* = 1, NN and NNN stand for nearest- 
neighbour and next-nearest-neighbour and each spin pair 
is counted once in the summation. The third term is 
an easy axis anisotropy, needed to select the ordering 
spin direction along c, and to generate a spin gap in the 
spectrum. In the absence of a clear physical picture of 
the origin of the anisotropy, we have chosen a single- 
ion form for this term. We anticipate that an exchange 
anisotropy [5JYSfS^ would have qualitatively similar 
effects. Both in-plane exchanges are antiferromagnetic 
Ji, J2 > whereas the inter-layer couplings are taken 
to be ferromagnetic Jz,Jz < 0, to stabihze the observed 
magnetic structure shown in Fig. 1151 

We have calculated the spin-wave dispersion and dy- 
namical correlations for the spin Hamiltonian ^ in the 
large S'-limit using the standard Holstein-Primakoff for- 
malism and details are given in Appcndix|X]with the dis- 
persion relations in eq. (|A2[) . To illustrate the effect of the 
various terms in the Hamiltonian we plot the dispersion 
for a triangular lattice with no anisotropy in Fig. 1 10b .) 
(dashed line). The linear spin- wave dispersion has zero 
modes not only at the magnetic Bragg peak position M 
for the specific ordered magnetic domain, but also at soft 
points such as M' in Fig. [S] which are Bragg peak po- 
sitions for the symmetry-equivalent domains (there are 
three possible magnetic domains related by a three-fold 
rotation around the c-axis). 

The presence of such unphysical gapless modes is a 
consequence of the macroscopic ground state degener- 
acy at the classical level and including 1/S corrections 
[Fig. [TUk ) solid line] generates a gap at the soft points 
but preserves a gapless Goldstone mode at the Bragg 
wavevectors.— The situation is completely different in the 
presence of an easy-axis anisotropy which generates a gap 
everywhere [see Fig.llOb)]. the gap at the Bragg wavevec- 
tor increases very rapidly as a power law, whereas at 
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FIG. 10; (Color online) Dispersion relation for the stacked triangular antiferromagnetic layers within LSWT (dashed lines) 
and SWT+l/S (solid lines) for various values of the parameters. The dispersion is plotted along high-symmetry directions in 
reciprocal space (dashed line path in Fig. |8^), top labels are high symmetry points in the Brillouin zone, a) Triangular lattice 
with antiferromagnetic couplings Ji = 1.32 meV, J2 = O.lSJi and no anisotropy. 1/ S quantum corrections have a large effect at 
the soft point M'. b) Same as a) but with easy-axis anisotropy I? = 1.78 meV, 1/S corrections are much smaller, c) Triangular 
layers with easy-axis anisotropy and coupled by NN interlayer exchange Jz = —0.3 meV. This leads to two incommensurate 
spin- wave modes with very little inter-layer dispersion near the gap minimum (last panel in the row), not consistent with data, 
d) Triangular layers with anisotropy coupled by second- layer exchange = —0.14 meV, this model gives the best fit to the 
data. Dashed (solid) line shows best fit using LSWT (SWT + l/S) with parameter values given in Table U 



the soft point increases only linearly with anisotropy, so 
above some threshold anisotropy value the minimum gap 
is no longer at the magnetic Bragg wavevector M, but 
moves to the soft point M'. 

Coupling the layers by a NN exchange < (ferro- 
magnetic) has the effect of splitting the dispersion into 
two modes [see Fig. [TUh)]. corresponding to acoustic and 
optic magnon modes between the two layers in the unit 
cell. The gap minimum is no longer at the commensu- 
rate soft point M' but at incommensurate positions sym- 
metrically displaced from M' along the magnetic stripe 
direction, see Fig. [TOb). The incommensurate position of 
the gap minimum varies linearly with \Jz\ for small 
and it originates physically from the fact that hopping 
of magnons is somewhat frustrated as an up spin in the 
bottom layer interacts through with two up and one 
down spin in the layers above and below, see Fig. 115b ). 
Moreover, there is very little interlayer dispersion of the 
gap minimum, i.e. between soft points M' with different 
L- values [see Fig.fTUb) last panel], also due to frustration 
in the magnon interlayer hopping at those wavevectors. 

Both of these features, incommensurate minima in the 
dispersion displaced (in-plane) away from M' and very 
weak inter-layer dispersion of the gap minimum are in- 
consistent with the data, which shows minima in the scat- 



tering at wavevectors very close to those corresponding 
to the commensurate M' point and a significant (~ 10% 
of the total bandwidth) interlayer dispersion of the gap 
minimum. This suggests that the observed interlayer dis- 
persion is not due to the nearest-neighbour interlayer ex- 
change Jz but is due to another inter-layer coupling. 

To find a possible explanation for the observed inter- 
layer dispersion we now assume = and consider 
the effects of a ferromagnetic second-layer coupling 
straight up along the c-axis, see Fig. [2b)- The resulting 
dispersion is plotted in Fig. [TUH) dashed hne, there is a 
single mode with essentially the same features as for de- 
coupled triangular layers [see Fig. [TUb)] with global min- 
ima at the commensurate soft point M' but now there is a 
modulation in the energy as a function of L- value to first 
order in . This model can account for all observed dis- 
persions in the boundaries of the powder-averaged spec- 
trum and the calculation for the best fit parameter set 
is plotted in Fig. [TTk-b) to be compared with data in 
Fig.HJi-c). 

In the calculation we used the dynamical correlations 
in eq. (jA3[) . included polarization factor and magnetic 
form factor as in eq. (|A5p . and performed a spherical 
average. An overall intensity scale factor C ~ 0.57 gives 
intensities comparable to data, we do not attach special 
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FIG. 11: (Color online) Powder averaged spin-wave spectrum for easy-axis triangular layers with/without (top/bottom) inter- 
layer couplings. Left panels to be compared with Fig. 3^) and right panels to Fig. The parameters are Ji = 1.32 meV, 
J2 = O.I5J1, T) — 1.78 meV, Jz = 0, and Ji = —0.14 meV (top panels) or = meV (bottom panels). The calculations 
contain an additive flat background. The various lines overplotted onto the intensity color maps show the dispersion along the 
directions from the origin F point to the closest soft points at M'(l/2,-l/2,0) (and the L — 1 and L — 2 versions of this) and 
M"(l,-l/2,0), with minima at wavevectors indicated by the vertical arrows in the right panels. The dispersion at the lower 
boundary of the scattering in b) is very similar to the data in Fig. |4j;) and are attributed to the weak interlayer couplings (no 
interlayer dispersion is present in d) calculated for decoupled layers). The dotted lines at low-Q mark the left edge of the region 
probed by the experiments. 



significance to this value except that it is comparable to 
the LSWT prediction of C = 1 . The intensity distribution 
is also reasonably well captured by the model as shown by 
the typical scans in Fig. [121 any differences compared to 
the data may be due to interactions beyond the minimal 
spin Hamiltonian considered. 

To obtain quantitative values for the model param- 
eters we have used the fact that various parts of the 
spectrum are mostly sensitive to only one parameter, 
in particular the gap depends only on anisotropy via 
Am' = 251? (1 — jg), the upper boundary of the spec- 
trum is mostly sensitive to the main coupling Ji, the 
magnitude of the low-energy dispersion between the gap 
minima with L = 1 and 2 is given by 45*17^1. The re- 
maining parameter J2 is extracted from the slope of the 
low-Q dispersion up to the first gap minimum. Specifi- 
cally, a fit to scan B at the minimum gap [Fig.H^B] gives 



the anisotropy V, scan C probes the inter-layer disper- 
sion so gives J^, fits to scans A and E upper boundary 
[Fig. n^E] give the main coupling Ji and fits to scan A 
lower boundary give the secondary coupling J2. 

The best fit values for the parameters obtained this 
way are shown in Table [D The calculated powder- 
averaged spectrum is shown in Fig. Illb -b) and compares 
well with the data in Fig. S^-c). To emphasize the sen- 
sitivity of the data to the 3D couplings we also plot in 
Fig. [TTb-d) the powder-averaged spectrum assuming de- 
coupled layers, which reproduces the bandwidth, gap and 
slope of the low-Q dispersion up to the first gap minimum 
but not the low-energy dispersion ('^ 10% of bandwidth) 
between subsequent minimum gap wavevectors. 

In the above analysis we neglected the NN interlayer 
coupling Jz since the powder data is not very sensitive 
to its presence, in particular a finite Jz produces only a 
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FIG. 12: (Color online) Comparison between data (filled points) and the linear spin-wave model for easy-axis triangular layers 
(solid lines) . Dashed lines show the estimated non-magnetic background and parameters are as in Fig. [TTk ) . 



TABLE I: Table of best fit values for easy-axis Heisenberg 
model applied to AgNi02. In both cases fits are obtained by 
comparing the spin wave dispersion calculated within linear 
spin wave theory, and including leading interaction correc- 
tions, with the envelope of the powder-averaged spin struc- 
ture factor S{Q, E) measured by inelastic neutron scattering. 
We set 5=1 and Jz = in all fits. All exchange parameters 
are given in meV. 





S 


Ji 


J2 


V 


J. 


JL 




LSWT 


1 


1.32(5) 


0.15(3) 


1.78(5) 





-0.14(2) 




SWT + l/S 


1 


1.35 


0.15 


1.59 





-0.14 





negligible inter-layer dispersion of the gap minimum (dis- 
persion appears only at order J^/I? for small Jz)- Using 
the fact that a finite Jz produces a small incommensurate 
shift of the minimum gap wavevector an upper bound for 
the maximum Jz that would still be consistent with the 
data can be estimated as | Jz|< 0.25 meV. 

Within the minimal model we considered in eq. (jl]) the 
estimated values of the interaction parameters would give 
a Curie- Weiss constant expected in a high-temperature 
susceptibility measurement on a powder sample as^ 



k^e = --S{S + l)(6Ji + 6J2 + 2J'z) = -5.9meV 
o 



This is smaller than, but comparable to, the experi- 
mental value of UbO = —9.2 meV extracted from high- 
temperature susceptibility data-^ and the difference may 
be due to other interactions beyond the minimal model 
we have considered to explain the dispersive boundaries 
of the powder data. 



C. Spin- wave theory with 1/5 corrections for 
easy-axis triangular layers 

In the absence of anisotropy the 2D triangular Ji — J2 
antiferromagnet in the coUinear phase has been shown to 
exhibit a strong renormalization of the dispersion rela- 
tion due to quantum fluctuations. Including the effects 
of fluctuations to order l/S in spin wave theory^ gener- 
ates a gap at the soft points where the LSWT calcula- 
tion predicts unphysical gapless modes, as illustrated in 
Fig.Iinii) sohd line is SWT-hl/S" and dashed is LSWT. 

In Appendix IB] we have calculated the l/S corrections 
to the dispersion relation for the triangular antiferromag- 
net when an easy-axis anisotropy is also present and the 
results are plotted in Fig. [TOb). The energy at the soft 
point is now only marginally renormalized up in energy 
because the easy-axis anisotropy generates a gap in the 
spectrum and this reduces significantly the zero-point 
quantum fluctuation effects. Using the 2D dispersion 
relation at order l/S and including also the inter-layer 
coupling the best parameter values obtained from com- 
parison to the data are very similar to those found using 
LSWT and are shown in Table [H The dispersion relation 
at order 1/5* for these parameters is plotted in Fig. [TUti) 
(solid line) and is very similar to the LSWT dispersion for 
the uncorrected parameters (dashed line). The powder- 
averaged spectrum is also similar to Figs. [TTk-b) for the 
uncorrected parameters and is therefore not reproduced 
here. 

In converging to the above set of parameter values 
for the microscopic Hamiltonian we have used two con- 
straints imposed by the data i) the inter-layer disper- 
sion is relatively small, ~ 10% of the bandwidth, as 
inferred from the observed small dispersion in the low- 
energy boundary of the magnetic inelastic scattering, and 
ii) the spin gap is relatively large, ^ 30% of the band- 
width, which in turn implies a relatively large easy-axis 
anisotropy. Specifically, subtracting the estimated inter- 
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FIG. 13: (Color online) Constraints on the parameters J2/J1 
and V/Ji determined from the spin gap measured in inelas- 
tic neutron scattering experiments. Gaps at M and M' are 
calculated within SWT+l/S' discussed in Appendix [B] and 
normalized to the maximum energy for spin excitations. The 
upper solid line divides solutions where the calculated gap at 
M' is larger than that seen in experiment, from those where 
it is smaller. The point J2/J1 ~ 0.15, V/Ji ~ 1.2 which 
gives the best fit overall fit to experiment is marked with an 
open circle. The region of small T> for which neither the gap 
at M nor M' are large enough to fit experiment is bounded 
by a dashed line. Linear spin wave calculations give similar 
lines, slightly shifted. [The small region of parameter space 
T>/Ji > 0, J2/J1 ^ 1/8 for which spin wave theory is ill- 
conditioned is shown by the red hatched area.] 



layer dispersion magnitude from the observed maximum 
energy leads us to search for a 2D dispersion shape where 
the gap energy divided by the maximum, i.e. the reduced 
gap A/wniax, is ~ 25%. 

The dispersion of the microscopic Hamiltonian in eq. 
(|4l) has the minimum either at the Bragg peak position 
M, or at the soft point M'. This places a strong constraint 
on model parameters, since both reduced gaps must be at 
or above the threshold value of 25%, with at least one gap 
precisely at the threshold. In Fig. [13] we map out these 
constraints for in ( J2/ Ji,!?/ Ji) parameter space, with 
gaps calculated using spin wave theory with 1/S correc- 
tions, as described in Appendix IB] Possible solutions lie 
on the solid line dividing solutions where the reduced gap 
AM'/'^max > 25% from those where AM'/^max < 25%. 

On this line other features in the data such as the ob- 
served slope of the low-Q dispersion further narrow down 
possible parameter values, leading to the best fit values. 
So within the minimal microscopic Hamiltonian consid- 
ered the requirement for quantitative agreement for the 
absolute gap relative to the maximum energy requires a 
relatively large anisotropy, and in this case the disper- 
sion minimum is at the soft point M' as illustrated in 
Fig. [Mi) (solid line). 
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FIG. 14: (Color online) Ground state phase diagram of the 
easy-axis Heisenberg model on a triangular lattice as a func- 
tion of J2/J1 and T>/Ji, calculated from mean field theory. 
The parameters J2/J1 « 0.15, V/Ji « 1.3 deduced from lin- 
ear spin wave theory fits to inelastic neutron scattering data 
place AgNi02 within the coUinear AF phase, not far from the 
boundary with a coplanar 3-sublattice state. 



V. DISCUSSION AND CONCLUSIONS 

We have reported inelastic neutron scattering measure- 
ments of the powder- aver aged spin dynamics in the lay- 
ered hexagonal antiferromagnet 2iJ-AgNi02, which has 
a coUinear alternating stripe ordered ground state. A 
broad band of magnetic inelastic scattering is observed 
above a finite gap indicating strongly dispersive magnetic 
excitations. The observed modulations in the boundaries 
of the magnetic scattering could be well explained by 
a gapped and predominantly two-dimensional dispersion 
relation, suggesting magnetically weakly-coupled trian- 
gular layers. For a quantitative analysis we have con- 
sidered a minimal spin Hamiltonian in eq. ^ for the 
localized Nil ions (5 = 1) located on an ideal triangular 
lattice and have neglected the coupling to the itinerant 
electrons on the surrounding Ni sites. Within this lo- 
calized model we considered first- and second-neighbour 
antiferromagnetic couplings Ji and J2 = 0.15Ji in the 
triangular layers, a strong easy axis anisotropy modelled 
by a single ion term V ~ 1.3Ji and weak interlayer cou- 
plings — — O.lJi and have found that this minimal 
model can account quantitatively for all the observed 
wavevector-dependence of the boundaries of the powder- 
averaged spectrum. 

The deduced values of the magnetic interactions pro- 
vide a natural explanation for the stability of the ob- 
served coUinear alternating stripe magnetic order. The 
large value for the easy-axis anisotropy strongly modifies 
the physics of an antiferromagnet on a triangular lat- 
tice. Firstly, it stabilizes coUinear order at the expense 
of coplanar states like the 120° spiral ground state of a 
simplest nearest-neighbour Heisenberg model. And sec- 
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ondly, it opens a gap to spin excitations which strongly 
suppresses quantum zero-point fluctuations. These are 
essential for the "order from disorder" selection of 
coUinear order in the isotropic J1-J2 model considered 
by Chubukov and Jolicocur'^. But in the easy- axis J1-J2 
model relevant to AgNi02, fluctuations lead only to small 
quantitative corrections to the spin wave dispersion. 

We have calculated these corrections explicitly to 
0{1/S) and find that, for most purposes, they can safely 
be neglected as they only lead to a very modest renormal- 
ization of the spin wave dispersion — and the parameters 
inferred from it — as illustrated in Table HI Competing 
interactions like J2 also play a secondary role, select- 
ing the observed form of coUinear order from the vast 
manifold of Ising ground states on a triangular lattice. 
The extent of the coUinear AF phase as a function of 
(J2/ Ji,!?/ Ji) can easily be estimated within mean field 
theory. The results of this analysis are shown in Fig. [TH 
The parameter set deduced above, with strong anisotropy 
V but relatively small second neighbour interaction J2 
places AgNi02 within the coUinear AF phase — as re- 
quired — but not far from the boundary with a coplanar 
3-sublattice state. 

At present, the origin of this strong easy-axis 
anisotropy — which may be of single-ion or exchange 
character — is unclear. AgNi02 is a highly anisotropic 
material, and the crystal fields at the magnetic sites will 
reflect this. However, relativistic (i.e. spin-orbit cou- 
pling) effects in Ni are not usually strong, so the absolute 
size of anisotropy is surprising. The more massive Ag 
ions, which have stronger spin-orbit coupling, are non- 
magnetic, and thought to play a negligible role in in- 
teractions between the Ni^+ sites. However, given the 
partially itinerant nature of Ni electrons, a simple local- 
moment model like eq. ^ is probably unlikely ever to 
offer a complete description of the magnetism in AgNi02. 

Further band-structure calculations based on the 
model developed in Ref. '5' could potentially test the pro- 
posed minimal Hamiltonian and give further insight into 
the microscopic origin of the exchange interactions and 
easy-axis anisotropy. Of particular interest would be the 
contribution of the itinerant Ni2 and Ni3 sites which sur- 
round the strongly-magnetic Nil sites to magnetic ex- 
citations, and their role in mediating interactions be- 
tween localized Nil sites. Even though the dispersion 
of the boundaries of the powder-averaged data are well 
explained by eq. there are discrepancies in the in- 
tensity distribution as a function of energy and wavevec- 
tor [see Fig. fl^E]. and these may contain information 
about the itinerant electrons. It is also interesting to 
ask whether the weak, but long-range interlayer interac- 
tion , could be of RKKY origin ? This is a promising 
avenue for future research, and recent transport exper- 
iments on AgNi02 in high magnetic field reinforce the 
idea that charge carriers couple strongly to the magnetic 
order.— 

Even without these complications, the minimal micro- 
scopic model we have proposed has a number of very 




FIG. 15: (Color online) Magnetic structure of 2ff-AgNi02. 
a) In each layer spins form alternating ferromagnetic stripes 
(dashed lines), ± symbols indicate the projection of the spin 
moments along the c-axis. Stripes are parallel between ad- 
jacent layers (thick black symbols for layer 1 and faint gray 
for layer 2) but have an offset because of the structural ar- 
rangement of Nil ions. Labels 1-4 refer to the four magnetic 
sublattices (two for each layer) in the magnetic unit cell (solid 
contour), which is doubled along a-axis compared to the crys- 
tallographic unit cell, b) 3D view of the magnetic structure 
showing the orientation of the spin moments (thick arrows) 
of the Nil ions (dark brown spheres) in the two layers in the 
unit cell. Light gray spheres are Ni2 and Ni3 ions, assumed 
non-magnetic. The box shows the 3D magnetic unit cell. 



interesting features. In particular it predicts a spin wave 
dispersion with global minima not at magnetic Bragg 
wavevectors but at symmetry-related soft points. This 
prediction could be tested by future inelastic neutron 
scattering experiments on single crystals or by ESR ex- 
periments, which probe the excitations at the F point, 
which by periodicity have the same energy as those at the 
magnetic Bragg wavevector. We note that the unusual 
dispersion relation we propose with global minima at 
multiple soft points may lead to non-trivial magnetically- 
ordered phases in applied magnetic fields that overcome 
the spin gap and future studies of the phase diagram 
in high applied magnetic field along the easy-axis could 
show a rich behaviour that would be fruitful to explore 
experimentally'^ and theoreticallyjiii 
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APPENDIX A: LINEAR SPIN- WAVE 
DISPERSIONS FOR COUPLED TRIANGULAR 
LAYERS 

Here we outline the derivation of the hnear spin-wave 
dispersion relations for a model of stacked antiferromag- 
netic triangular layers with collinear stripe order. The 
experimentally-determined magnetic structure in 2H- 
AgNi02 is shown in Fig. [T21 In each triangular Ni layer 
the ordered spins form alternating ferromagnetic rows 
described by two magnetic sublattices (up/down), lead- 
ing to a total of 4 magnetic sublattices for the two lay- 
ers in the magnetic unit cell (box in Fig. [T5b). After a 
Holstein-Primakoff transformation from spin to magnon 
operators and Fourier transformation the spin Hamilto- 
nian in eq. ([4]) becomes 



where Eq is a constant and terms higher than quadratic 
are neglected. The sum extends over all wavevectors q 
in the first magnetic Brillouin zone. The q dependence 
of the operator matrix X and of the Hamiltonian matrix 
H is implicit. The operator matrix is given by 



Xt=[c 



, /3- 



where a, /3, e and 7 refer to the magnetic sublattices 1, 2, 
3, 4 (in this order) and creates a plane-wave magnon 
mode on sublattice 1, and so on. The Hamiltonian matrix 
in this operator basis is 



H = 



A B C D* 
B* A D C 
C D* A B 
D C B* A 



(Al) 



where 



A =2S [Ji cos(27rfc) + Ji + J2 cos(27r(2ft. + k)) + J2 

+J', cos(27r0 - - J,] + 2SV (l - ^) 
B = ASJz cos(7rfc) cos(7rZ) C 
C =45JiCos(7rA;)cos(7r(2/i + fc)) 
cos(37rfc) cos(7r(2/i + k)) 
2SJ, cos(7r0 

-{2h+k)Tri/3 



D 

c 



where (h, fc, I) are components of the wavevector q in 
units of the reciprocal lattice of the hexagonal unit cell 
shown in Fig.[2]3). Note that the very last term in func- 
tion A above, 25*1? x ^ formally comes in spin-wave the- 
ory as a higher order {1/S) correction term, but we have 
included it here following common convention)^ since it 
originates from a quadratic term in boson operators. 

Diagonalizing the Hamiltonian in eq. (|A1[) using stan- 
dard methods gives two doubly-degenerate modes with 



dispersions given by 

(w±)2 ^ A^ + BB* -C^ - DD*± 



sgn(B/C) x/4\AB - CD*\^ - \B*D* - BD\^ 

(A2) 

where ± stands for optic and acoustic modes between the 
two layers and sgn(a:) is the sign function. The above 
expressions for the dispersion relations also hold by peri- 
odicity for a general wavevector Q = {h, k, I) outside the 
first magnetic Brillouin zone [the two modes are degen- 
erate when B — and the prefactor sgn(_B/C) in front 
of the square root is required to ensure continuity of the 
two distinct modes on both sides of a degenerate point] . 

For decoupled triangular layers with no anisotropy 
{Jz = J'z — T> = the dispersion is plotted in Fig. 
[TUk ) (dashed line) and for the general case is plotted 
in Fig. [TOhXdashed and dash-dotted lines). The finite 
easy-axis anisotropy 2? > leads to a gap in the spec- 
trum. The gap above the magnetic Bragg peaks in- 
creases as a power law with increasing anisotropy o^j^ = 

25^2? (4 Ji + 4J2 + 2Jz + , where V ^ V{l-^). 

However, the gap is near-linear in anisotropy at the soft 

V — 2J,] , so in the limit of 



point M', Wjyr, 



25* 



^2 



weak couplings between the layers the minimum gap is 
always at the soft point M' and not at the Bragg wavevec- 
tor M [see Fig.[T35)]. 

The dynamical correlations (per Nil spin) are obtained 

as 



S^^ {Q,uj)=C§ 



\W{-LJ+)+X{-u;+)+Y{~uj+) + Z{-uj+)\ 



\W{-L^-)+X{-UJ-) + Y{-UJ-) + Z{-UJ-)\'^ 



\W(uj^)+XiLO-)+Y(uj-} + ZiLO.)\^ 



G{lo~u^) (A3) 



where G{liJ — u!±) is a Gaussian with a finite width to 
model the instrumental resolution, we use the shorthand 
notation ijj± = uIq and the functions W , X, Y, Z and N 
are given by 



W{lo) = -{A + lu){A^ + BB* -C^ 
2ABB* - C{B*D* + BD) 



DD* -uj^)+ 



X{uj) = C{A^ + BB* -C^ + DD* - w^)- 
A{B*D* + BD) - uj{B*D* - BD) 

Y{uj) = B* [{A + uj)^ - BB* + 6*2] - 
2C{A + u;)D + BD^ 

Z{lu) ^ D{A^ ~ DD* -uj^) + B*^D*- 
2AB*C 
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WW* 



XX* - YY* 



ZZ* 



The overall intensity prefactor C in eq. (|A3p is 1 in LSWT 
and is considered here as a variable parameter in the com- 
parison with the experimental data in order to account 
(in a first approximation) for a possible overall intensity 
renormalization compared to LSWT. We note that the 
above expressions simplify considerably in the absence of 
the interlayer coupling J^, in that case the two magnon 
modes in eq. (IA2[) become degenerate with dispersion 
UJQ = — and the dynamical correlations in eq. 
(|A3|) become 

S^--{Q,u;)=C^^—^G{u;-loq). (A4) 

For completeness we also quote the obtained transfor- 
mation matrix S to the normal operator basis X' = S^^X 
in which the Hamiltonian is diagonal 



S = 



W{uj-) W{lu+) Wi-LU^) Wi-LU+) 

F(c._) Y{uj+) n-uj^) Y{-uj+) 

X{lj^) X{lj+) X{-uj+) 

Z{lO^) Z[tjJ-^) Z{ — LU^) Z{ — LU^) 



where we use the shorthand notation W{uj) — 
W{uj)/ y/ N{u!) and so on. This matrix satisfies the eigen- 
value equation gHS = SgH' and the normalization con- 
dition SgS^ = g, where g is the commutator matrix for 
the operator basis defined by 



x{x*y - (x*x^ ) 



10 
10 
0-10 
-1 



and 
















































H' = 



is the Hamiltonian representation in the normal operator 
basis X'. 

The one-magnon neutron scattering intensity including 
the polarization factor and the magnetic form factor is 



5(g,^) = (7rof (^1 + glj (|/(Q)) S'^^{Q,u) 

(A5) 

where {•^To) =290.6 mbarns/sr is a conversion factor to 
bring the intensity into absolute units of mbarns/meV/sr, 
f{Q) is the magnetic form factor for Ni^+ ions and we 
take the (7-factor g = 2. Qz is the component of the 
wavevector transfer Q along the c-axis. A numerical av- 
erage of eq. (jA5l) over a spherical distribution of orienta- 
tions of the wavevector transfer Q was performed in order 
to obtain the orientational-averaged intensity S{Q,lu) to 
be compared with the measured powder data. 



APPENDIX B: SPIN- WAVE THEORY WITH 1/S 
QUANTUM CORRECTIONS FOR THE 
EASY- AXIS Ji - J2 TRIANGULAR 
ANTIFERROMAGNET 

Here we outline calculations of corrections to spin wave 
dispersion arising from interactions between magnons at 
0(1/5). Calculations to this order have previously been 
carried out by Chubukov and Jolicoeur for the isotropic 
J1-J2 Heisenberg model on a triangular lattice^. Here 
we extend this work to the case of a finite easy-axis 
anisotropy V. To make contact with this earlier work 
we use here an orthogonal coordinate system with x 
along the magnetic stripes and y transverse to stripes 
in plane as in Fig. [5)d) and label the 2D wavevectors as 
) — {Qxa,Qya) where Qx.y are orthogonal compo- 
nents (in A~^) and a is the hexagonal lattice parameter. 

In the usual linear, semi-classical approximation (i.e. 
neglecting terms which contain more than two Bose op- 
erators) we obtain a dispersion 



with coefficients 



Ak = 25 Ji ^1 + cos kx + a + a cos Vsfc, 
Bk = 45 Ji cos ^ (cos ^ -f a cos ^ ) (Bl) 



251? 



where a = J2/J1 and V ^ V{1 - 1/(25)). We note 
that our units differ by an overall factor of two from 
those in Ref. H, and that single-ion easy-axis anisotropy 
contributes terms to the linearized Hamiltonian which 
are formally of order 0(1/5) relative to usual LSWT. 

The leading effect of interactions between spin waves, 
treated at a one-loop (i.e. Harteee-Fock) level, is to renor- 
malize the coefficients of this dispersion to give 



efc = ViAk + Mfc)2 - {Bk + 5Bkf 



where 



5Ak = 
SBk = 



E 



2Ji 



A. 



B. 



2L0r, 



Bk^-^Gik,p) [B2) 



2uj. 



Here, N is the total number of spins in the triangular 
plane, and the sum on internal momenta p extends over 
the entire nuclear Brillouin zone (large thick line hexagon 
in Fig. Hi). Following Ref. we write Z = Z/{2SJi) 
(Z = A, B, Lu) and introduce the functions 



F{k,p) (1 - cosfcj;)(l - cospj,) - 2(1 + ( 
+a{l — cosViky){l — cos Vipy) 



2V 

X 



16 



G{k,p) — cos 



V^ky + y/5py + Px 



y^ky kx 's/Spy Px 

■ cos COS 



+a 



y/Sky + 3fc^ \/3py + Spx 
cos 'hi COS - 



2 2 
V3fcy -3kx V^Py - 



(B3) 



Anisotropy D enters SA^ explicitly through F{k,p), 
but its main effect is to introduce a gap into the non- 
interacting dispersion (D^, which acts as an infra-red cut- 
off in all denominators. In order not to mix different 
orders in 1/S, we set V ^ V in calculations of SAk and 
5Bk. 

To obtain quantitative results for comparison with 
experiment, the sums in eq. (|B2[) were evaluated nu- 
merically. The resulting corrections to the linear spin- 
wave dispersion are illustrated first for the case of no 
anisotropy in Fig. [TOb ) (solid line). In this case we ex- 
actly reproduce the results of Ref. : LSWT (dashed 
line) predicts a gapless mode at the soft point M' (as an 
indirect consequence of the macroscopic classical ground 
state degeneracy), however quantum corrections lift the 
classical degeneracy by stabilizing the coUinear stripe or- 
der and generating a spontaneous gap at the soft points 
while maintaining a gapless Goldstone mode at the mag- 
netic Bragg wavevectors. 

In the presence of anisotropy the situation is some- 
what different. Firstly the coUinear order observed in 
AgNi02 is now stable at a classical level for a wide range 
of I? and J2 (c.f. Fig. [M]). And secondly, spin excita- 
tions are now gapped, even in the absence of fluctuation 
effects [see Fig. [TOb)]. Expanding about colHnear order, 
at a semiclassical level, the gap at the Bragg wavevectors 
M = (0,27r/\/3) grows as 



Am«25^I?(4Ji+4J2+P) 

whereas at the soft points M' = (tt, tt/Vs) it increases 
only linearly as 2S'D. As a result, for finite V, the min- 
imum in the spin wave dispersion occurs on the zone 
boundary at the soft point M', and not at the magnetic 
ordering vector M. Once 1/5* corrections are included, a 
gap is generated dynamically at M' for T> — 0, and there 
is a threshold anisotropy value "Dc above which the global 
minimum gap moves from M to the soft point M'. 

For parameters relevant to 2i7-AgNi02 (see Table |l| 
the gap to spin excitations is large throughout the BZ, 
with a minimum at M' {Vc « 0.21 meV for this pa- 
rameter set). Quantum fluctuations lead only to small 



quantitative corrections to the dispersion — compare the 
dashed to the solid line in Fig.ITOb) — and fits to LSWT 
are therefore justified a posteriori. A further indication 
that quantum fluctuations are small is the fact that the 
zero-point reduction of the sublattice magnetization 



N ^ 2u>„ 



p 



(B4) 



is calculated to be small : 5ms = 0.038 for the parameters 
given in Table |T1 

We note that this is not true in the limit of 
small anisotropy 1) ^ and small NNN couplings 
J2/J1 1/8, in which case the system is close to a quan- 
tum critical point. In this case zero-point fluctuations are 
very strong and the correction to the sublattice magne- 
tization diverges within LSWT. A consequence of this is 
that l/S corrections drive the gap 



Ar ~ i^yiJiS + J2S)V[S - 1/2 ~ Sms] (B5) 

at the r point k = (0, 0) imaginary for a small but fi- 
nite region of parameters V/Ji > 0, J2/J1 > 1/8 — 
c.f. Fig. [131 Accurate treatment of this small V regime 
would require a more sophisticated self-consistent spin 
wave theory, which is not warranted for AgNi02. 

For the purposes of comparing the spin-wave calcula- 
tions at 0(1/5*) with the envelope of scattering in (|Q|, uj) 
it is sufficient to replace A^ , and lu^ by their corrected 
values in the expression for the scattering intensity 



_ ^SAk 
~ ^2 



5{iv-ujk). (B6) 



Finally, the interlayer coupling can easily be in- 
cluded in these calculations by modifying Eq. (|B1[) as 
follows 

Ak Afc + 25*4 (cos fc^-1). 

Here kz — QzC where Qz is the c-axis wavevector com- 
ponent in A~^. At order 0{\/ S) the interlayer couplings 
modify Eq. (|B3|) to read 



F{k,p) ^ F{k,p) + ^il-cos2kz){l~cos2pz), 

however since this change is relatively very small for sim- 
plicity we neglect it in the analysis and only consider 
the explicit effects of through changes in the dis- 
persion relation at linear order in SWT. The resulting 
three-dimensional spin wave dispersion and corrections 
are shown in Fig. [TUti). 
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